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Satellite-based measurements of the Earth's atmosphere and surface are very important 
because they help us understand our planet's climate, monitor global air quality, and predict the 
weather. Almost all of these measurements are affected by clouds. Some instruments are 
designed specifically to study how clouds impact climate. For other measurements, clouds can 
either be a nuisance or they may actually help us to extract information about gases in the 
atmosphere. In all cases, it is important to understand exactly how clouds impact the satellite 
observations. 

Ozone is an important constituent of the Earth's atmosphere, and it is a focus of several 
space-based instruments. It acts as a protective shield by absorbing ultraviolet rays high in the 
atmosphere. But ozone in the atmosphere near the Earth’s surface can also be harmful to life. It 
damages lung tissue when inhaled and can create visible scars on plants. It is important to be able 
to determine how much ozone is in the upper atmosphere where it is crucial to our survival and 
how much is in the lower atmosphere where it is considered to be a pollutant. 

Satellites are extremely useful for measuring ozone globally. However, satellite instruments 
do not directly sample the Earth's atmosphere. Instead, they make measurements in different 
wavelengths of light either reflected from the sun by the atmosphere, clouds, and surface or 
emitted as heat. The measured wavelengths include colors that we can see, invisible light that 
can bum our skin, and heat (including microwaves) from the atmosphere, surface, and clouds. 
Because clouds are good reflectors of light, they can shield the lower part of the atmosphere 
from satellite instruments. We can use this property and the fact that clouds vary in height to 
slice up the atmosphere and tell us where exactly the ozone is. But first we must understand 
precisely how clouds affect the incoming sunlight. 

There are currently 5 satellites flying in a formation; They observe the same regions of the 
Earth's atmosphere within minutes of each other. This formation is known as the A-train because 
tire first satellite is named Aqua and the caboose is called Aura. Both Aqua and Aura are part of 
N ASA’s Earth Observing System. One of the middle cars, called Parasol, carries an instrument 
that can determine the height of a cloud using the absorption of sunlight by atmospheric oxygen. 
Aura has an instrument that can make similar measurements using two completely independent 
techniques. This paper shows that all three techniques provide similar estimates of the cloud 
height. Some of the small differences can be traced to features of the individual retrieval 
algorithms. This comparison serves as a means of validating our algorithms. 
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Abstract. The cloud pressures determined by three' different algorithms, 
operating on reflectances measured by two space-borne instruments in the 
“A” train, are compared with each other. The retrieval algorithms are based 
on absorption in the oxygen A-band near 760 nm, absorption by a collision 
induced absorption in oxygen near 477 nm, and the filling in of Fraunhofer 
lines by rotational Raman scattering. The first algorithm operates on data 
collected by the POLDER instrument on board PARASOL, while the lat- 
ter two operate on data from the OMI instrument on board Aura. The satel- 
lites sample the same air mass within about 15 minutes. 

Using one month of data, the cloud pressures from the three algorithms 
are found to show a similar behavior, with correlation coefficients larger than 
0.85 between the data sets for thick clouds. The average differences in the 
cloud pressure are also small, between 2 and 45 hPa, for the whole data set. 
For optically thin to medium thick clouds, the cloud pressure the distribu- 
tion found by POLDER is very similar to that found by OMI using the O 2 -O 2 
absorption. Somewhat larger differences are found for very thick clouds, and 
we hypothesise that the strong absorption in the oxygen A-band causes the 
POLDER instrument to retrieve lower pressures for those scenes. 
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1. Introduction 

Clouds have a large influence on the transfer of radiation in the atmosphere. This makes 
clouds important in climate studies and for trace gas retrievals in passive remote sensing. 
For climate studies several properties are needed: particle phase, particle radius, cloud 
liquid- or ice-water content, cloud optical thickness, and cloud (top) pressure or cloud 
(top) temperature. These are usually observed using a combination of wavelength bands 
in the visible and thermal infra-red part of the spectrum. For the cloud correction of trace 
gas retrievals from UV /VIS reflectance spectra two much simpler cloud parameters are 
commonly used: an effective cloud fraction c eff and a cloud pressure p c . These parameters 
are found from a fit of the observed top-of-atmosphere reflectance, and the strength of a 
height-sensitive spectral feature. In the present article we compare cloud pressure data 
from two satellite instruments flying in the “A” train, using one month of data with global 
coverage. 

This comparison includes three cloud products: cloud pressure derived from the O 2 A- 
band absorption at 760 nm, cloud pressure derived from O 2 -O 2 absorption at 477 nm and 
cloud pressure derived from the filling in of Fraunhofer lines by rotational Raman scat- 
tering at 350 nm. The first is observed by the POLDER (Polarization and Directionality 
of the Earth’s Reflectances) instrument on PARASOL (Polarization and Anisotropy of 
Reflectances for Atmospheric Sciences coupled with Observations from a Lidar), the lat- 
ter two are observed from OMI (Ozone Monitoring Instrument) on Aura. The POLDER 
instrument is specifically designed to study cloud and aerosol properties from space, while 
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21 OMI is designed to measure high resolution reflectance spectra to perform atmospheric 

22 composition measurements. 

23 The structure of this paper is as follows. The next section briefly described the two 

24 instruments, followed by a section on the cloud retrieval algorithms. Next is a short section 

25 on matching measurements from OMI to measurements from PARASOL, followed by a 

26 description of the actual comparison results. We end with a discussion of the similarities 

27 and differences we observe, and a brief discussion of future improvements. 

2. Description of the instruments 

2 s Both Aura and PARASOL are part of the so called “A” train, a series of satellites 

29 carrying Earth observation instruments. Near the front of the train is the PARASOL 

30 satellite with its POLDER instrument, which will be described in brief detail in section 2.1. 

3 1 The last satellite in the A train is Aura, which carries four instruments, including OMI. 

32 This instrument is briefly described in section 2.2. Both instruments sample the same 

33 part of the atmosphere within approximately 15 minutes. PARASOL has a local equator 

34 crossing time of about 13:30, Aura crosses the equator at about 13:45. 

2.1. Description of PARASOL/POLDER instrument 

35 PARASOL is flying in formation with Aqua and Aura (NASA), CALIPSO 

36 (NASA/CNES) and CloudSat (NASA/CSA) as part of the A train. The PARASOL scien- 

37 tific objectives are to characterize the radiative and microphysical properties of clouds and 
sa aerosols using as best as possible the data complementarities from the different sensors 

39 on board the A train. PARASOL is carrying a wide-field imaging radiometer/polarimeter 

40 called POLDER. POLDER is designed to measure the directionality and polarization 
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4 1 of light reflected by the Earth-atmosphere system. The POLDER instrument is exten- 

42 sively described by Deschamps et al. [1994], It is a digital camera with a two-dimensional 

43 (274 x 242 pixels) charged coupled device (CCD) detector array, wide field of view tele- 

44 centric optics and a rotating wheel carrying spectral and polarized filters (see Fig. 1). 

45 Similar POLDER instruments have already flown aboard the Japanese ADEOS-1 (1996— 

46 1997) and ADEOS-2 (2003) platforms. Contrary to those first versions of POLDER, for 

47 the PARASOL version the telecentric optics array has been turned 90 degrees to favor 

48 multidirectional viewing over daily global coverage. When the satellite passes over a tar- 

49 get, up to 16 observations are realized (up to 14 with the previous configuration). The 
so swath is now 1600 km (across track) corresponding to a maximum field of view of 114°. 

51 A 490 nm polarized channel was also put in place of the 443 nm one. Moreover a 1020 nm 

52 waveband has been added to conduct observations for comparison with data acquired 

53 by the lidar on CALIPSO. The spectral bands and the central wavelengths of POLDER 

54 aboard PARASOL are reported in Table 1. 

55 This instrument presents original features since it is not only multispectral but also mul- 
se tidirectional and multipolarization. Algorithms dedicated to “Earth Radiation Budget, 
57 Water Vapor, and Clouds” were developed, taking into account these capabilities [Buriez 
sa et al., 1997]. More particularly, the multi-polarization capability allows determining the 
59 cloud thermodynamic phase and the cloud top pressure, the multi-directionality improves 
so the derivation of the cloud optical thickness and the estimate of the reflected flux, whereas 
si the multi-spectrality allows deriving the cloud middle pressure and the clear-sky water 

62 vapor content. Daily products and monthly syntheses are produced at 20 km resolution 

63 (after cloud detection performed at full resolution, 6 km, and for every direction). The 
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64 data archive starts from March 4 th , 2005, and PARASOL is still operational at present 
es time. 

2.2. Description of OMI on Aura 

66 The Ozone Monitoring Instrument (OMI) is a contribution of the Netherlands’ Agency 
6 ? for Aerospace Programs (NIVR) in collaboration with the Finnish Meteorological Institute 
es (FMI) to NASA’s EOS Aura mission. OMI will continue the TOMS satellite data record 

69 for total ozone and other atmospheric parameters related to ozone chemistry and climate. 

70 The OMI instrument employs hyperspectral imaging in a pushbroom mode to observe 

71 solar backscattered radiation in the visible and ultraviolet. The observed spectra cover 

72 the wavelength range 270 nm to 500 nm, with a spectral resolution of 0.42- 0.63 nm. The 

73 swath is wide enough to allow for global coverage in one day (14 orbits), with a spatial 

74 resolution of 13 x 24 km 2 for nadir observations. The spectral range and resolution of 

75 OMI allows for the retrieval of column amounts of atmospheric trace gases, like O 3 , NO 2 , 

76 SO 2 , BrO, HCHO, cloud detection is needed to correct those trace gas retrievals for the 
„ presence of clouds. 

78 OMI uses two 2-dimensional charged coupled device (CCD) detector arrays, one for the 

79 UY wavelength range (270- 350 nm) and the second one for visible wavelengths (350- 

80 500 nm). On either CCD, one dimension is used for the separate wavelengths, while the 
ai perpendicular dimension is used for the 60 across track positions (see Fig. 2). Unlike 

82 GOME, Sciamachy and GOME-2, OMI has no scanning mirror and its response is made 

83 independent of the polarization of the detected radiation with the use of a polarization 

84 scrambler. A detailed description of the OMI instrument and its science objectives can 

85 be found in Levelt et al. [2006a, b]. 
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3. Short overview of the cloud height retrieval algorithms 
a« Two of the retrieval algorithms use absorption of radiation by oxygen to determine the 
a? height of clouds in the atmosphere, while the third uses the amount of rotational Raman 
aa scattering observed from the filling in of the Fraunhofer lines in the solar spectrum to 

89 determine the cloud pressure. They all use reflected sunlight, rather than thermal infra- 

90 red emissions from clouds, as is done in most meteorological satellite retrieval techniques 
si for cloud top temperature and cloud top pressure. The oxygen absorption feature used in 

92 the first two algorithms is rather different, as is the spectral resolution of both instruments. 

3.1. POLDER cloud pressure retrieval using the oxygen A-band at 760 nm 

93 Two different methods were developed to retrieve cloud pressure from POLDER data. 

94 The first one (cloud Rayleigh pressure) is based on the analysis of polarized reflected 

95 light at 490 nm, and is not discussed further in the present article. The second one 
ss (cloud oxygen pressure) uses the ratio of the two POLDER radiances measured in the 

97 oxygen A-band near 763 nm [ Buriez et al., 1997] . Cloud oxygen pressure po 2 is determined 

98 from differential absorption between the radiances measured in the channels centered 

99 at 763 nm (narrow band) and 765 nm (wide band) respectively (see Fig. 3). The 

100 and Ry 65 radiances are first corrected for gaseous absorption of ozone and water vapor, 

101 then the measured oxygen transmittance Tq 2 is obtained from the ratio of f ?763 and 
loa f? 765 - All the gaseous transmissions are derived from simulations using a line-by-line 

103 model [Scott, 1974], The spectroscopic database used for the absorption cross sections is 

10 4 HITRAN 2004 [ Rothman et al., 2005]. In the first step, the influence of the surface albedo 
los is neglected. An apparent pressure p app is inferred by assuming that the atmosphere 
we behaves as a pure absorbing medium overlying a perfect cloud reflector located at pressure 
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10 ? p app . In practice, p app is calculated from a polynomial function of Tq 2 and the geometric 
wo air-mass factor M = 1/ cos 6 + 1/ cos/9 0 - The coefficients of the polynomials are fitted 
m from line-by-line calculations. 

no Because of enhanced oxygen absorption due to the effects of surface reflection and 
in multiple scattering inside the cloud, the apparent pressure p app is almost always higher 
in than the cloud top pressure. For example, even for optically thick clouds, large differences 
in (typically 200 hPa) were observed between POLDER-1 apparent pressures and cloud top 
in pressures derived from the brightness temperatures measured in the 11 pm channel of 
ns METEOSAT [Vanbauce et al, 1998]. Comparable differences were observed between 
us the apparent pressure and the Rayleigh pressure derived from POLDER polarization 
11 ? measurements [Parol et al., 1999]. The apparent pressure can even be higher than the 
ns cloud base pressure when a great amount of photons reaches the surface before being 
ns reflected back to space, that is in the case of a thin cloud layer above a bright surface. 

120 Cloud oxygen pressure po 2 is determined from the apparent pressure by removing the 

1 2 1 surface contribution. This correction is only realized for pixels over land surface, because 

122 the ocean reflectance is low at 765 nm and therefore the surface influence is negligible. 

123 Over sea-surface only viewing directions outside the sun-glint are retained. The scheme 

124 of the cloud oxygen pressure algorithm is given in Fig. 4. The starting point is that the 

125 oxygen A-band corresponds to strong absorption lines for which the oxygen transmission 
m Tq 2 can be treated by means of a random band model [Goody, 1964]: 


Tq 2 = exp(-C'v / Mp ap p) 


( 1 ) 
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127 where M is the geometric air mass factor and C a constant depending on spectroscopic 

128 data. Considering that this transmission can be decomposed in a term corresponding to 

129 the light directly reflected by the cloud and a term corresponding to the light reflected 
no after reaching the surface, the surface-corrected oxygen pressure can be written after some 
131 approximations (see Vanbauce et al. [2003] for details) in: 


Papp “t (u 1 )Psurface 

Po 2 = 

r 


( 2 ) 


132 where r is the fraction of photons directly reflected by the cloud and p surface is the surface 

133 pressure. The fraction of photons reflected by the cloud, r, is calculated using r — 

134 R® 65 /R 765 where R 7 q 5 is the reflectance measured by POLDER at 765 nm after correction 
us for gaseous absorption and Rj 65 is the reflectance that would be measured if in addition 
us the surface was black. p s „ rface is obtained from the ECMWF (European Center for Medium 

137 range Weather Forecasts) analysis. In the operational algorithm, po 2 is calculated only 

138 for cloudy pixels with optical thickness larger than 3.5. 

139 From comparisons of P OLDER- 1 cloud oxygen pressure and ARM/MMCR [Clothiaux 
mo et al. , 2000] cloud boundaries pressures, po 2 appears to indicate the cloud middle pressure 
141 rather than the cloud top pressure [ Vanbauce et al. , 2003] . 


3.2. OMI cloud pressure retrieval using the collision induced absorption at 
477 nm 

w2 Only a brief overview of the OMI O2-O2 cloud model and cloud retrieval algorithm 

143 will be given here, since they are described in considerable detail in Sneep et al. [2007b] 

144 and Acarreta et al. [2004], All atmospheric oxygen absorption bands (A, B, and 7 bands, 
us the oxygen transition a 1 A g (u = i ) <— X 3 T,~(v = 0) for i = 0,1,2, respectively) fall 
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165 


outside the wavelength range of OMI. This means that the FRESCO method for cloud 
height detection [ Koelemeijer et al., 2001], which is used for GOME and Sciamachy is 
not readily available for OMI. However, oxygen has several collision induced absorption 
(CIA) features within the OMI wavelength range, and they may be used instead. In these 
CIA features two oxygen molecules jointly absorb a single photon, and each fly away 
from the collision in an (electronically) excited state. The strongest of these CIA features 
within the OMI wavelength range is found at 477 nm, see for instance Greenblatt et al. 
[1990]. Because the absorption cross section of O 2 -O 2 scales with the squared number 
density of oxygen, rather than directly with the oxygen number density as is the case 
for the oxygen A-band, some care is needed to correctly retrieve a cloud pressure from 
observations at 477 nm, and some different biases may be expected, compared to FRESCO 
or the POLDER oxygen cloud pressure. 

A DOAS (Differential Optical Absorption Spectroscopy [Platt, 1994]) fit of the OMI re- 
flectance spectrum between 460 and 490 nm is used to determine the slant column amount 
of O 2 -O 2 . This value, combined with the viewing- and solar geometry and surface condi- 
tions, is used to find the cloud pressure with the aid of a lookup table. The lookup table 
was produced with the DAK (Doubling Adding KNMI [de Haan et al., 1987; Stammes, 
2001]) radiative transfer model, using a Lambertian surface with albedo 0.8 as the cloud 
model. Simulations have shown that the pressure of the cloud retrieved by this method is 
at about the mid-level of the cloud [Sneep et al, 2007b], even for optically thick clouds. 

3.3. OMI cloud pressure retrieval using the filling in of Fbaunhofer lines by 
rotational Raman scattering at 350 nm 
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Rotational-Raman scattering (RRS) causes filling-in and depletion of solar Fraunhofer 
lines throughout the ultraviolet in the observed backscattered Earth radiance (normalized 
by the solar irradiance) [e.g. Joiner et al, 1995]. This property was first used to retrieve 
an effective cloud pressure by Joiner and Bhartia [1995], Spectral fitting methods that 
exploit the high-frequency spectral structure of RRS have been applied to hyperspectral 
instruments such as GOME and OMI [Joiner et al., 2004; Vasilkov et al., 2004; Joiner 
and Vassilkov, 2006]. The latter reference contains a description of a soft-calibration 
procedure that is used to remove scan position-dependent biases (i.e. striping) from the 
retrieved cloud pressures. 

The OMI RRS algorithm is currently implemented with the same cloud model as the 
OMI O 2 -O 2 cloud retrieval algorithm, as described in section 3.4. There are two sets of 
products based on separate sets of assumptions applied to this model: The first set of 
products is included for historical reasons using a cloud albedo of 0.4 that produces an 
effective cloud fraction close to the MODIS geometrical cloud fraction. A second set is 
produced assuming a cloud albedo of 0.8 that gives cloud pressures closer to the physical 
cloud top at the lower cloud fractions. The latter set of products (called ‘CloudPressure- 
for03’ and ‘CloudFractionfor03’ in the OMCLDRR product files) is the one that will be 
used throughout this paper. 

These products are generated assuming a fixed surface albedo of 0.15 that was chosen 
to be consistent with the OMI total ozone retrieval based on the Total Ozone Mapping 
Spectrometer (TOMS) version 8 algorithm. This value is known to be higher than the 
actual surface albedo under most conditions but was designed to account for aerosol and 
small amounts of low-level cloud in the OMI TOMS-V8. In an off-line study, we have 
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applied the assumption of a 0.05 surface albedo to the OMCLDRR algorithm. We found 
that this assumption brings the cloud pressures into closer agreement with the OMI O 2 -O 2 
cloud algorithm especially at the lower cloud fractions. 

3.4. Differences in the cloud models used by POLDER and OMI 

Both OMI cloud products use basically the same cloud model, which is the same as 
the cloud model used in FRESCO [Koelemeijer et al, 2001]. The cloud is represented 
by a Lambertian surface with albedo 0.8, no light is transmitted through the cloud. The 
scene is partially covered by the model cloud with an effective cloud fraction c eff , so that 
the top-of-atmosphere reflectance agrees with the observed reflectance. The albedo of the 
model cloud is so high that most scenes have an effective cloud fraction less than one; the 
missing transmission of this model cloud is compensated by the large cloud-free part of 
the pixel. Comparisons with simulations of scattering clouds have shown that the albedo 
of 0.8 is a suitable value for this model cloud [Koelemeijer and Stammes , 1999; Wang 
et al., 2006; Vasilkov et al . , 2007]. The cloud pressure is adjusted so that the retrieved 
cloud shows the same amount of signal (either O 2 -O 2 slant column, or amount of Ring 
effect) as the observation. 

The POLDER cloud model is different from the OMI cloud model, namely a scattering 
and transmitting cloud. Here the retrieval is limited to cloudy subpixels (6 x 6 km), where 
there is complete cloud cover with an optical thickness of 3.5 or larger. Over sea, where 
the surface is very dark at 760 nm, the cloud optical thickness is used as a threshold value 
in determining the cloud pressure. Over land, where the surface can be very bright at 
760 nm, especially over vegetation, the cloud optical thickness is used both for selection 
and correction of p app . The cloud pressures measured from different viewing angles are 
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averaged, and then the results for the cloudy sub-pixels are combined with a cloud cover . 
weighted mean into the final cloud pressure at 18 x 18 km 2 pixels. 

4. Matching individual scenes in OMI and PARASOL 

The pixels on which POLDER reports the cloud pressure are 18 x 18 km 2 , comparable 
to the OMI nadir pixel size of 13 x 24km 2 . For this reason a one-to-one mapping between 
the two datasets was chosen, with a single PARASOL scene compared to one OMI scene. 
The PARASOL data is stored on a non-rectangular grid, and functions exist to map a 
(latitude, longitude) coordinate pair onto this grid. For each OMI pixel the matching 
PARASOL pixel is looked up, and stored on the OMI grid for later comparison. For this 
article a special dataset was prepared where each orbit is stored in a separate file, rather 
than the standard single day in an orbit. This was done to avoid overlap of successive 
orbits at higher latitudes. 

5. Comparison results 

For this comparison a total of 383 orbits were used (OMI orbit numbers 9986 to 10422, 
PARASOL repeat cycle 34, orbit 219 to cycle 36, orbit 189), covering most of June 2006. 
The two instruments sample the same part of the atmosphere within about 15 minutes. 
The measurements were filtered to exclude pixels over a bright surface by excluding snow 
or ice covered surfaces. For these scenes it is known that the contrast between cloud cover 
and the surface is too low to properly distinguish clouds from the background, leading to 
an incorrect effective cloud fraction [Sneep et al., 2007b], and therefore an ill-determined 
cloud pressure. Furthermore, the data was filtered to exclude pixels with a POLDER 
cloud cover less than 95 %, and pixels where the rotational Raman effective cloud fraction 
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is less than 0.2, because the rotational Raman algorithm switches to a different cloud 
model in those cases. The OMI rotational Raman scattering cloud product comes in two 
flavors; here the “cloud pressure for O3” was used exclusively. 

Histograms showing the global distribution of cloud pressures from the three retrieval 
methods are shown in Fig. 5 separately for scenes over land and sea. Over sea a bi-modal 
pressure distribution is found, while over land only a single mode is observed. Although 
the overall shape of the distribution of cloud pressures is very similar, some differences 
can be seen. To investigate where these differences occur, separate histograms are made 
for small (0.2 < c eff < 0.4) and large (c eti > 0.8) effective cloud fractions (from the OMI 
O2--O2 algorithm) , shown in Fig. 6. The distributions of the differences between the three 
cloud pressures are shown in Fig. 7. These observations will be discussed in section 6. 

Scatter plots of all combinations of the three parameters are shown in Fig. 8, again 
separated for land and sea. The correlation coefficient p and the slope from a straight 
line fit including the errors in both data sets, following Press et al. [2003, section 15.3], 
are listed in each of the sub-figures. 

Fig. 9 shows the correlation coefficients, the median difference, and the 66 % quantile 
width between all three data sets over land and over sea as a function of the effective 
cloud fraction. An increase in correlation with increasing c eS is seen for land and sea. The 
median difference shows some interesting behaviour which will be dicussed in section 6. 
The results are summarized in table 2. 

6. Discussion 

The three cloud pressure products are in good to excellent agreement, with average 
differences between them that are well within the stated accuracy of those products. 
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253 From other comparisons and model studies [Vanbauce et al, 1998; Koelemeijer et al, 

254 2001; Vanbauce et al. , 2003; Sneep et al, 2007b; Vasilkov et al. , 2007] it was already clear 

255 that the cloud pressure derived from visible or near infrared reflectance spectra is well 

256 within the cloud, and probably close to the mid-pressure level. This is in stark contrast to 

257 thermal infrared observations, where the cloud top pressure is retrieved. An exception to 

258 this rule is the cloud Rayleigh pressure from POLDER, where the degree of polarization at 

259 490 nm is used, and the underlying assumption is that a cloud will scramble all polarization 

260 signal, yielding the top of the cloud layer, sometimes even above the cloud top pressure 
S61 found by a thermal infrared instrument like MODIS [ Parol et al, 2006]. 

262 Not only are the average differences small, the correlation between the data sets is high 

263 and the slope observed in the scatter plots is reasonably close to 1 , giving confidence in 

264 all algorithms involved. With measurements that are in such good agreement, there are 

265 details that tend to stand out, and those details will be discussed below. 

266 From the distributions shown in Fig. 5, in particular over sea, one could conclude that 

267 the OMI O 2 -O 2 cloud pressure retrieval is less sensitive for low pressure clouds than the 

268 0 2 A-band retrieval from PARASOL. One might expect that this is caused by the pressure 

269 dependence of the absorption strength of the collision induced absorption (<7o 2 -o 2 °c V 2 ) • 

270 On the other hand, the rotational Raman scattering product does not have a similar 

271 pressure dependence, and yet it shows a similar behavior at low pressures compared to 

272 the OMI O 2 -O 2 cloud pressures. Model studies presented in Sneep et al. [2007b] indicated 

273 that the expected influence of the quadratic pressure dependence of the absorption cross 

274 section is limited to approximately 40 hPa, which can not explain the median difference 
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275 of ^lOOhPa found here for thick clouds. Because the differences are most clearly seen 

276 over sea, we limited the next few steps to that subset. 

277 Inspection of Fig. 5 for pixels over sea shows that for clouds at low pressures the PARA- 

278 SOL O2 A-band algorithm retrieves smaller pressures than the OMI O2-O2 and RRS algo- 

279 rithms. A similar effect can be seen in Fig. 6 for pixels with a large effective cloud fraction. 

280 In these cases we deal presumably with convective clouds with the cloud top located at 

281 low pressures. The OMI RRS and O2 0 2 algorithms need to put the Lambertian cloud 

282 at relatively high pressures, corresponding to pressures deep inside the scattering cloud, 

283 to reproduce the measured signal [Vasilkov et al, 2007 ]. In contrast, the O2 A-band algo- 

284 rithm can put the perfect reflector at lower pressures, closer to the cloud top, to reproduce 

285 the measured signal. Due to the relatively strong absorption in the O2 A-band photons in 

286 this band may not penetrate as deeply inside the scattering cloud, while photons in the 

287 weakly absorbing 0 2 A 3 2 band and photons affected by Raman scattering penetrate deep 

288 inside the scattering cloud. Therefore, the 0 2 -C >2 and RRS algorithms retrieve higher 

289 pressures than the O2 A-band algorithm for these clouds. For optically thin clouds, which 

290 are probably also geometrically thin, photons can penetrate the entire cloud for all of 

29 1 the three algorithms. Therefore, similar distributions are found for the O2 A-band and 

292 the O2-O2 band for small effective cloud fractions in Fig. 6. The deviating behaviour of 

293 RRS for thin clouds is believed to be caused by the assumed value of the surface albedo. 

294 In Sneep et al [ 2007 a] it is shown that the cloud pressures retrieved by the RRS method 

295 are much closer to the O2-O2 cloud pressures when an improved surface albedo is used 

296 for the RRS method. 
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297 From a qualitative comparison with CloudSat radar profiles, we hypothesise that the 

298 more frequent occurrence of clouds between 700 and 750 hPa in RRS, seen most clearly 

299 in the thick cloud distribution shown in Fig. 6, is caused by a combination of effects: 1) 
am the surface albedo assumption in RRS, which causes it to be too low, 2 ) effects of the 
sox cloud model used, which could well be different for both OMI cloud products since there 

302 is more Rayleigh scattering at the wavelengths used for RRS, and differences in the way 

303 multi-layer cloud decks are handled. The presence of sun glint has opposing effects on 
sew both OMI products, causing a shoft towards low pressures for RRS and a shift towards the 

305 surface for O2-O2. The effect of sun glint on the present analysis was investigated, and 

306 while the correlation between the two OMI cloud pressures improved slightly at low cloud 

307 fractions, no significant changes in the statistical results were obeserved. More research, 

308 including radiative transfer calculations in geometrically thick clouds and multiple cloud 

309 decks, are needed to understand the differences between the algorithms. 

7 . Conclusions and outlook 

310 The cloud pressures retrieved from OMI and POLDER measurements using oxygen ab- 
3u sorption or the amount of rotational Raman scattering to determine the cloud height find 

312 remarkably similar cloud heights. In general the cloud pressure measured by these meth- 

313 ods is much higher than the cloud pressure derived from thermal infrared measurements. 

314 Model studies and comparisons with ground based radar profiles [ Vanbauce et al, 1998 ; 

315 Koelemeijer et al, 2001 ; Vanbauce et al., 2003 ; Sneep et al., 2007 b; Vasilkov et al, 2007 ] 

316 suggest that the cloud pressures retrieved here indicate the mid-level of the cloud layer. 

317 Despite the good agreement, there are some differences visible between the three al- 

318 gorithms, due to different sensitivities and different assumptions used at various stages 
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319 in the retrieval. The OMI O2-O2 algorithm uses a monthly surface albedo climatology 

320 derived from GOME measurements at 1 ° x 1 . 25 °, while the rotational Raman scattering 

321 algorithm uses a fixed value for the surface albedo of 0.15 which comes from the TOMS 

322 heritage. In a future version both will switch to a surface albedo climatology derived from 

323 OMI measurements at 0 . 25 ° x 0 . 25 °. This will affect the cloud fraction most directly, 

324 but a change in effective cloud fraction will change the cloud pressure because the same 

325 strength of the spectral feature needs to be explained. 

326 The strength of the oxygen A-band leads to a different sensitivity to the cloud optical- 

327 and geometrical thickness when compared to the much weaker oxygen collision induced 
32 s absorption at 477 nm or rotational Raman scattering near 350 nm. This difference affects 

329 the retrieved cloud pressure for scenes with a high effective cloud fraction, where POLDER 

33 0 retrieves a pressure closer to the cloud top than the other two algorithms. 
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CCD detector 



Figure 1. The measurement principle of POLDER on PARASOL. 



Figure 2. The measurement principle of OMI. 

Table 1. The spectral bands in POLDER on PARASOL. Channels labeled with (P) measure 


polarization. 


Channel 

Bandwidth 

Rationale 

443 nm 

20 nm 

Ocean color applications 

490 nm (P) 

20 nm 

Cloud properties, Aerosol retrieval 

565 nm 

20 nm 

Calipso lidar at 532 nm 

670 nm (P) 

20 nm 

Aerosol retrieval, Cloud properties 

763 nm 

10 nm 

Cloud oxygen pressure by differential 

765 nm 

40 nm 

absorption in oxygen A-band 

865 nm (P) 

40 nm 

Aerosol retrieval, Cloud properties 

910 nm 

20 nm 

Water vapor retrieval 

1020 nm 

20 nm 

Calipso lidar at 1064 nm, Aerosol retrieval 
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Figure 3. POLDER/PARASOL filter transmissions in the narrow and wide bands centered at 
763 nm and 765 nm, respectively, together with atmospheric transmission in the oxygen A-band 
region. 



Figure 4. Scheme of the POLDER cloud oxygen pressure algorithm. 
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1000 800 600 400 200 

p c for c eff > 0.2 /hPa 


Figure 5. The distributions of cloud pressures from the OMI O 2 -O 2 , the OMI rotational 
Raman scattering, and the POLDER on PARASOL O 2 A-band products, for scenes over land 
(top) and sea (bottom). 




Figure 6. The distribution. of cloud pressures from the OMI O 2 -O 2 , the OMI rotational Raman 
scattering, and the POLDER on PARASOL O 2 A-band products, over sea for scenes with a large 
effective cloud fraction (top) and scenes with a small effective cloud fraction (bottom) . 


DRAFT 


March 19, 2007, 4:29pm 


DRAFT 


SNEEP ET AL.: OXYGEN CLOUD PRESSURES FROM OMI AND POLDER 


X- 27 



Figure 7. The distribution of differences in the cloud pressure between the O2-O2 cloud 
pressure, the rotational Raman scattering, both from OMI on EOS Aura and the oxygen cloud 
pressure from POLDER on PARASOL for colocated scenes over sea, for scenes with a large 
effective cloud fraction (top) and scenes with a small effective cloud fraction (bottom). 
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Figure 8. Probability distribution of the cloud pressure determined from OMI and PARASOL. 
The contours represent the densest area in the- scatter plot, with the contours containing 10%, 
30%, 60%, 90%, and 99% of all points, going to progressively lighter colors, for each of the 
three combinations of two algorithms. The data is shown separately for land and sea surfaces. 
The dotted line in each of the plots are the x = y relation, the drawn line is the result of an 
orthogonal regression analysis, the slope of which is printed in each plot. 
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Figure 9. Correlation, 66% central quantile width and median difference between all three 


combinations of cloud pressure products, over both land (drawn lines) and sea (dashed lines), 


plotted as a function of the effective cloud fraction. The measurements were grouped by c off , 


from 0.2 to 0.4, from 0.4 to 0.6, from 0.6 to 0.8, and 0.8 and larger. 
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Table 2. Some statistical parameters describing the differences of the co-located cloud pressure 
retrievals. The difference is the product listed at the top minus the product listed at the left, 
the slope is for the product listed at the top projected on the horizontal axis. This is for pixels 
over land and sea combined, filtered to include only pixels with c off > 0.5. 
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Ap c = 45 hPa 
(j(Ap c ) = 74hPa 
p = 0.93 
slope = 1.19 

A p c = 2 hPa 
a(Ap c ) = 93hPa 

p = 0.88 

slope = 1.32 

OMI 

O 2 -O 2 

A p c = —45 hPa 
cr(Ap c ) = 74hPa 
p = 0.93 
slope = 0.84 


Ap c = — 44hPa 
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p = 0.92 
slope = 1.09 
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p = 0.88 

slope = 0.76 
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slope = 0.92 
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